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The determination of the quantum properties of a single 
mode radiation field by heterodyne or double homodyne de- 
tection is studied. The realistic case of not fully efficient pho- 
todetectors is considered. It is shown that a large amount of 
quite precise information is avalaible whereas the complete- 
ness of such information is also discussed. Some examples 
are given and the special case of states expressed as a finite 
superposition of number states is considered in some detail. 
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I. INTRODUCTION 



developments have extended the method to data coming 
from inefficient detectors These various procedures, 
unitarily referred to as quantum tomography, provide a 
nice and powerful tools for investigating the quantum 
properties of radiation field leading to the maximum in- 
formation achievable on the measured state. However, 
they require the detection of many field quadratures x^. , 
j =,!,..., N and thus a lot of repeated measurements 
on the state under examination. A systematic approach 
to precision of quantum tomography is in progress 
however, a preliminary study ||7[ has indicated that its 
determination of some field properties can be very noisy 
relative to the direct measurements of the same quanti- 
ties. 
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In order to get information about a quantum state one 
has to measure some observable. A question immedi- 
ately arises: is this information complete ? Namely, does 
it fully specify the quantum state under examination ? 
The answer is generally negative: the measurement of 
only one observable show up only an aspect of a physical 
system whereas a complete description requires a deeper 
inspection. However, the measurement of several differ- 
ent observables could require considerable efforts. There- 
fore it is a matter of interest to compare the simplicity 
and the feasibility of a measurement, or a set of mea- 
surements, with the provided amount of information. In 
addition, one has to weigh up the precision of such an 
information. 

The complete description of a quantum state is con- 
tained in the density operator p, or for pure state in 
the wave function. Therefore, a measurement leading 
to the density matrix in some representation provides, in 
principle, all the knowable information about the mea- 
sured state. It has been shown theoretically [|l[ that 
the Wigner function Wo{a,d) of a field mode can 
be reconstructed starting from the homodyne measure- 
ments of a continuos set ip S [— tt, tt) of field quadratures 
Xip — l/2{ae^'^'^ + a^e*"^). Later this procedure has been 
applied to coherent and squeezed states Q using a fi- 
nite set of phases (pi, i = 1, ...27, upon a smoothing on 
experimental data. More recently a procedure suitable 
to finite sampling has been suggested for recovering ma- 
trix elements in the Fock representation 0] and latest 



II. REALISTIC HETERODYNE DETECTION 

In this paper a different approach to quantum state 
measurement will be investigated. The crucial remark 
is that the density matrix in the coherent state repre- 
sentation (a|p|Q;) is a positive definite function and thus 
can be directly measured for any quantum state of ra- 
diation field. Indeed, it is known that the so called 
Husimi Q-function VF_i(a, a) = (a|/5|o;) represents the 
outcomes probability distribution [|8|-pi| of both the het- 
erodyne ll^ and the double homodyne jlj] detection 
scheme when equipped with ideal photodetectors. On 
the other hand its smoothed versions 

Ws{a,a) ^ W.,{f3,P) exp{-2^^} , (1) 

emerge from realistic devices in which not fully efficient 
detectors are involved. The parameter s depends on 
quantum efficiency as |l3| , |l^ 



(2) 
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Starting from heterodyne, or equivalently from double 
homodyne [|l5|, detection we can evaluate some quantity 
O of interest as an average over the experimental distri- 
bution 

f (fa 

0={0)= J ^Ws{a,a) Ts [O] («, a) , (3) 

where J^s[0](a, a) is a non operatorial function related 
to the s-ordering, in the boson operator expansion, of the 
corresponding observable O |l^. Denoting by {0}^ the 
s-ordered form of an operator we have for example 
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a^a = {a^a}s + ^{s — 1) 
at2a2 ^ {at2a2}^ + 2{s - l){a^a}, + i(.s - if 

and thus 

Ts[n\{a, a) = jap + i(s - 1) 



(2s-l)|a|2 + is(s-l) 



(4) 



(5) 



for the mean photon number and for the number fluctu- 
ations. Similarly, we can investigate the squeezing prop- 
erties of the examined state by means of the s-ordering 
of the field quadrature fluctuations 



(a,a) = i[a2e-2»^+a2e2,^ 



2\a\ 



(6) 



and also checking the uncertainty product. This proce- 
dure can be generalized in order to evaluate any field 
correlation which is of interest. In fact, any t-ordered 
moment {a^"a"+''}f, with arbitrary t, can be written in 
terms of a finite number of s-ordered moment by means 
of the formula IIm 



^ {d + n)\ 



n — k 



(7) 



where also s is arbitrary. The expectation value of any t- 
ordered moment can thus be evaluated by an average over 
the different ordered distribution Ws{a,a). The special 
case in which the parameters t and s are choosen to be 
t = 1 and s = 1 — 2?7~^ is of interest as it allows to ob- 
tain the normal ordered field correlations (a^"a"+'') in 
terms of a finite numbers of heterodyne measured mo- 
ments {a^^a^^'^),-!- We have 



{d + n)\ ( n 



X / — l^i_2^-i(a,a) a" 
or in a more compact form 



k+d^k 



(8) 
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where L'^{x) denotes Laguerre polynomials. An interest- 
ing application of Eq. (^) lies in the evaluation of the 
expectation value (e'""^) of the shift operator. Starting 
from normal ordering 
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we have, in fact 



(10) 



d^a 



W^i-2,7-i(a,a) 



fc=0 



1 - e*"^ 
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Then, using properties of Laguerre polynomials [18|, we 
arrive at the formula (valid for < < arccos(l -^rp/2)] 



?/ - 1 + e 



i0 



d^a 



Wi^2r,-^ ("' ") 



X exp{ ^ } 



(12) 



Eq. (||) is also suitable for a reliable estimation of er- 
rors in the determination of the various expectation val- 
ues. In practical situation, in fact, one has at disposal 
a finite sample of Ws{a, a) and thus the integral in for- 
mula (^) has to replaced by the corresponding statistical 
sampling 



J G data 



(13) 



According to the law of large numbers O approaches (O) 
in the limit of infinite number of sampling data, whereas 
for finite sample we have a confidence interval O ± SO, 
with SO given by 



y J G data 



(14) 



Some examples of the present reconstruction procedure 
can be given by means of numerical simulations of realis- 
tic heterodyne detection. In Fig. |^ I report the results for 
the mean photon number (n) determination at different 
values of the quantum efficiency for coherent states of dif- 
ferent amplitudes. Fig. nk shows the results from hetero- 
dyne detection and Fig.nip that ones from a direct pho- 
todetection. The two determinations are also compared 
in Fig. |l|c. In making such a comparison the same num- 
ber of repeated measurements on the field have to be con- 
sidered. In a scheme of N repeated measurements of the 
quantity x, the accuracy Sx rescales as 5x oc iV~^/^. The 
proportionality costant generally depends on the kind of 
detection. If the outcomes x are distributed around the 
true value x according to the probability p(x\x), the er- 
ror for N repeated measurements is always bounded by 
the Cramer-Rao inequality [|l§ (5x > {FN)~'^/'^, F be- 
ing the Fisher information F — J dx [dxp{x\x)]'^ / p{x\x) . 
For Gaussian distributed data one has F — 1/cr^, with 
the variance of the distribution, and the lower bound 
for precision is achieved. In practical situations, in order 
to evaluate the precision Sx, one takes advantage of the 
central limit theorem [ p^ , which assures that the par- 
tial averages over a block of Nj, data is always Gaussian 
distributed around the global average over many blocks. 
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Thus, one evaluates precision by dividing the ensemble 
of data into subensembles, and then calculates the r.m.s. 
deviation of subensemble averages with respect to the 
global one. 

From Fig. ^ it is apparent that the method works also 
for low efficiency of the photodetectors even though the 
results are slightly more noisy than ones from direct de- 
tection. However, this level of introduced noise can be 
considered as admissible in sight of the further informa- 
tion available from the same data sample. Moreover, it 
has to be noticed (see Rcf. ^) that tomographic deter- 
mination of (h) introduces a very large amount of noise, 
even for unit quantum efficiency. Fig. |^ illustrates the 
application in recovering field fluctuations on a squeezed 
state and a number state for different values of the quan- 
tum efficiency. 

About the determination of the phase some further 
considerations are in order. Heterodyne detection, in 
fact, is by itself a phase detectors as the marginal dis- 
tribution 



(15) 



represents a phase distribution of the measured state 
|^ , pO| . We have thus at disposal not only the mean value 
of the phase and the fluctuations but also the whole prob- 
ability distribution. The distribution in Eq. ( |l5| ) does 
not coincide (even for 77 = 1 ) with the canonical phase 
distribution 



P(0) 
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i(n — m) 



Pn.r, 
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and it is generally broadened relative to (|l^). In Fig. ^ 
the phase distribution obtained for a squeezed vacuum is 
reported for experiments carried out with different val- 
ues of the quantum efficiency. The distributions broaden 
when the quantum efficiency decreases but the crucial 
information about phase bifurcation |^l| is not lost also 
for for very inefficient detectors. 

The results obtained until now can be summarized in 
the following assertions: i) starting from heterodyne de- 
tection many properties of the measured state can be 
determined at one go; ii) this determination is slightly 
more noisy than the corresponding ones from direct mea- 
surements of the same quantities one at times, even for 
unit quantum efficiency of the photodetectors. However, 
this behaviour is not unexpected as heterodyne detec- 
tion involves the joint measurement of non commuting 
observables, and thus the unavoidably addition of noise 
by first principles [ p2|j2^ . This is not the case of quan- 
tum tomography where each homodyne measurement is 
independently performed and noise is introduced by data 
processing itself. 



III. DENSITY MATRIX RECONSTRUCTION 

Let us now deal with the completeness of the informa- 
tion coming from heterodyne detection. Is it possible, as 
an example, to determine the whole number distribution 
? The matrix elements Pn+k,n in the Fock representation 
are the expectation values of the generalized projectors 



Pn,n+k = \n){n + k\ 



|0)(0| 



-.n+d 



y/{n + dy. ' 

Using the Louisell expansion of the vacuum jl^ 

e^l- ^ pi 



(17) 
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we can write Pn.n+k in terms of normal ordered moments 
1 



Pn,n+k = I ./ lini V — \ 

n + fc)! e^i" „ p! 



gt"+PQ^n+P+fc 



(19) 



Eq. ( p^ ) is suitable to apply Eq. (^. After a straighfor- 
ward calculation we get the reconstruction formula (|^) 
for the matrix elements 



Pn+k,n — 



id 



(n ^ k)\ ^-^ V ?7 

' q—n ^ ' 



X , — W^_^^-^{a,a)a^Ll{JM')■ (20) 



Unfortunately, Eq. (|20|) is not suitable for statistical sam- 
pling as the interchange of integration and summation 
is not mathematically allowed p4| . The analytical ex- 
pression of Ws(a, a) is needed and thus some smoothing 
procedure on sampled data is required, unavoidably in- 
troducing some a priori hypothesis on the state under 
examination ||2^. However, Eq. (|2^) is far from being 
a purely formal tool. In many situations, in fact, one 
has some information about the considered state and 
thus some suggestions on parametryzing Wigner func- 
tions. Moreover, the distributions Ws{a,Oi) for s < — 1 
are smoothed functions by themselves p6[ and generally 
do not exhibit large or fast oscillations. Therefore we may 
expect the smoothing not leading to a dramatic lost of 
information and, at the same time, to not introduce fake 
information. The completeness of information coming 
from heterodyne detection cannot, anyhow, be claimed 
in a general way. 

The reconstruction of the entire density matrix (in the 
Fock representation) and thus a complete description of 
the state is possible for the special case of states with a 
finite number of moments different from zero. Examples 
of such a states are provided by finite superpositions of 
number states 
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TV 



(21) 



n=0 



The latter can be produced in different manner in a high- 
Q cavities p6| , p7| and also by a special non linear inter- 
action lESi. If the moments a^"a'" vanish for n or m 



beyond a certain value the series in Eq. (19) are actually 
truncated and Eq. ( pO| ) becomes suitable to a statistical 
sampling 



Pn+k,n „ 



(-)" / n! [(fa 



Wi_2rt-^{a,a) a'' 




(22) 



p=0 



The value of N has to be choosen large enough to ensure 
the cancellation of any moment a^^^^ a'^^^ i,j = 0, 1, .... 
In practice one can start with a large value of N and then 
optimizing it by means of some stability criterion. In any 
case the precise value of N is not needed by the algo- 
rithm. In Table || 1 report the results of this procedure 
when applied to the state 



^±{\Q)+^\2)) 
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using photodetectors with an overall quantum efficiency 
equal to r/ = 0.9. The reliability of the method is appar- 
ent. The same degree of precision can be obtained with 
lower efficiency using a larger sample. 

The problem of reconstructing the density matrix of 
field-states with finite numbers occupation has been con- 
sidered also by other authors, in particular in the context 
of the so-called Pauli's phase retrieval problem, where 
two experimentally determined probability distributions 
are needed. Orlowsky and Paul proposed inj 2^ an algo- 
rithm to recover the entire wavefunction (|l), assuming 
that the moduli of the wavefunction are known in the 
position and momentum representation. Their method 
involves solving blocks of nonlinear equations after a de- 
composition of the wavefunction into Hermite polynomi- 
als. The resulting phase retrieval scheme is transparent, 
however it is recursive from above, namely it determines 
the coefficients 4'n from the highest index N to the low- 
est. In addition, the value N of the nonzero components 
of the wavefunction has to be known in advance. On 
the contrary, Eq. (^) allows recovering of the matrix 
elements Pn,m one by one as an average over the experi- 
mental distribution and the value of N is not needed by 
the algorithm. Bialynicka-Birula and Bialynicki-Birula 
in ||30|] considered the reconstruction problem starting 
from the knowledge of the photon number and phase 
(Pegg-Barnett) distributions. They suggested two dif- 
ferent algorithms based on recursive iterations of Fast 
Fourier Transform from the number representation to the 
phase domain. Their algorithms work only for pure states 



whereas the present one can also be applied in the general 
case. In fact, the only requirement for writing Eq. ( |20| ) 
in the sampling-suited form ( p2| ) is that the measured 
state possesses only a finite number of moments different 
from zero. This condition can obviously be fulfilled also 
by finite mixtures. It has also to be mentioned that a de- 
tection scheme for the Pegg-Barnett phase distribution 
has not been devised yet. Thus the phase distribution 
needed by the algorithms in | ]30| can only be inferred by 
other state measurement schemes such as homodyning or 
quantum tomography. 

Apart from the above considerations 1 consider the re- 
liability of the present method in evaluating the confi- 
dence interval on matrix elements determinations as its 
main advantage. 



IV. CONCLUSION 

In conclusion, quantum state measurement by hetero- 
dyne or double homodyne detectors has been shown to 
provide a large amount of quite precise information. It 
cannot lead to a complete specification of the measured 
state due to the singularity in some s-ordering (s < — 1) 
of operators. To the knowledge of the author it repre- 
sents, at current time, the best compromise between the 
conflicting necessity of precise and complete state mea- 
surement. 
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FIG. 1. Mean photon number determination by a simu- 
lated heterodyne and direct detections for different values of 
the quantum efficiency rj. In (a) the results for three different 
coherent states of amplitude a = 1 (circle), a — 2 (square) 
and a = 3 (triangle) are reported for heterodyne detections. 
In (b) are reported the results from direct detections with the 
same number of repeated measurements on the same coherent 
states. The confidence intervals of both the determinations 
are evaluated as usual, by dividing the whole sample of 10^ 
data in subensembles and then calculating r.m.s. deviation 
with respect to the global average (see text). In (c) the ac- 
curacy of the two kinds of determination is compared. The 
noise (in dB) added by heterodyne detection is shown as the 
ratio between the corresponding confidence intervals. 



FIG. 2. Simulated determination of field fiuctuations Ax^ 
for ip = 0,7r/2. Results for a squeezed state of total en- 
ergy (h) = 1 equally distributed between coherent amplitude 
and squeezing and a number state p = |1){1| are reported. 
The upper and the lower curves are for the squeezed state 
{if = 0,7r/2 respectively), the central one the result for num- 
ber state (isotropic field distribution). Samples of 10^ data 
are used and the confidence intervals are evaluated as in Fig. 





FIG. 3. Phase distribution from heterodyne detection for 
a squeezed vacuum with squeezing photons (n) — 1. The 
distributions are obtained with a sample of 10^ data for 
77 = 0.25,0.5,0.75,1.0. Broader distributions correspond to 
lower values of rj. 
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TABLE I. Reconstructed density matrix along with the confidence intervals for the state \ip) = l/\/2(|0) + i|2)). The 
experiment has been simulated for quantum efficiency ri = 0.9 using a sample of 10® data. 

.502 ±.024 (.004 - i.003) ± (.022 + i.021) (.001 - i493) ± (.037 + i038) 

(.004 + i.002) ± (.021 +i.022) -.003 ±.053 (-.003 ± i.002) ± (.018 ± z.018) 

(.001 ± i.493) ± (.037 - i.038) (-.003 - i.002) ± (.018 ± i.018) .500 ± .031 

-.004 ±.063 
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